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FOREWORD 

This  publication  was  prepared  under  contract 
by  the  UNITED  STATES  JOINT  PUBLICATIONS  RE¬ 
SEARCH  SERVICE,  a  federal  government  organi¬ 
zation  established  to  service  the  translation- 
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and  research  needs  of  the  various  government 
departments. 
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HEAT  WAVE  RADIATING  ENERGY  FROM  A  FRONT 

■  .  ;  .  -  ussr  -  ;  ;  .  '/ 

[Following  is  a  translation  of  an  article  by  E.  I.  Andri- 

anWri  in  Zhurnal  TekhniolieskojrFisiki  (Journal  of  Tech¬ 
nical  PhysicsTTTol.  1X1X7 No,  ll7”ioscow--Leningr-£d. 

November  1959*  pages  1363-1372* 1 

The  article  investigates  the  propagation  of  a  nonprcgressing 
heat  wave  radiating  energy  from  a  front#  It  examines  the  case  when 
the  travel  of  radiation  in  a  cold  gas  is  large  for  all  frequencies 
below  a  certain  critical  frequency  6C  and  small  for  higher  fre¬ 
quencies.  The  travel  of  quanta  in  the  heated  region  is  assumed  to. 
be  much  shorter  than  the  radius  of  the  wave  front,  so  that  the  radi¬ 
ant  transfer  of  energy  occurs  by  thermal  conductivity, 

When  a  heat  wave  of  high  temperature  propagates  through  a 
gas  its  molecules  split  Up  into  atoms,  the  latter  ionizing  to  form 
plasma  with  a  high  radiant  heat  conductivity  (l).  It  can  sometimes 
be  considered  that  the  radiation  travel  In  a  plasma  is  much  shorter 
than  the  radius  of  the  wave  front. 

The  coefficient  of  radiation  absorption  in  a  cold  gas  greatly 
depends  on  the  frequency  of  quanta  and  rapidly  approaches  zero  for 
long  wavelengths .  Therefore,  as  Kompaiieyets  showed,  it  can  be  con¬ 
sidered  that  only  those  quanta  leave  the  eave  front  whose  frequency 
is  lower  than  a  certain  Value  60  #.  The  problem  of  a  progressing 
heat  wave  propagating  from  a  point  in  which  there  is  a  sudden  liber¬ 
ation  of  a  quantity  of  energy  Q,  assuming  the  internal  energy  and 
the  heat  transfer  coefficient  to  be  exponentially  dependent  on  tem¬ 
perature  but  without  talcing  into  account  the  radiation  of  energy  from 
the  front,  was  investigated  by  ZelMovich  and  Kompaiieyets  (2). 

Let  us  examine  by  the  approximate  method  a  similar  problem, 
taking  into  account  the  radiation  of  energy  from  the  wave  front  and  , 
applying  the  relationships  of  moments  which  were  studied  in  detail 
by  Bareriblatt  (3)  in  the  theory  of  the  filtration  of  a  gas  in  a 
porous  medium,  and  by  Loytsyanskiy  (k)  in  the  theory  of  a  boundary 
layer.  :  ••••'  '''-r  :  v  ■-  •' 
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We  shall  use  for  our  study  a  particular  example*  when  the 
internal  energy  of  the  gas  is  proportional  to  the  temperature,  and 
the  length  of  the  free  travel  of  radiation  is  a  function  of  temper¬ 
ature  in  accordance  with  the  exponential  law  1  s  IqT^.  Let  us  note* 
however,  that  this  same  method  may  he  used  to  study  the  problem. when 
there  is  an  arbitrary  dependence  of  the . internal  energy  and  radiation 
travel  on  temperature.  0  ,  "  T-.  " '/• 

We  shall  consider  that  'the  cold  gas  passes  ail  frequencies 
lower  than  w  #  and  detains  frequencies  higher  than  (j*.  We  shall 
disregard  the  shielding  effect  on  radiation  due  to  the  heauing  of  the 
prefrontal  zone.  The  condition  for  energy,  balance  at  the  wave  front 
can  then  be  written  thus  . 


*  C° 


-  -  S  (T^), 


M. 


hw^d 


to 


hw 


k  =  q  /  U> 


CQ  Villas  , 
3k 


(1) 


|^2C2  e^’  -  li 


h  z  6.62i>*10“27  erg/sec j  k  =  1*38 ‘lO”3^  erg/degree;  a  —  heat  capa¬ 
city  per  unit  volume;  c  —  velocity  of  light*  the  index 
signifies  that  values  refer,  to  the  wave  front.  In  a  number  of  eases, 
S  can  approximate  the  exponential  depend  slice.  If  hw  <4~kT  tp,?  then; 
resolving  the  exponent  into  the  series  (Rayleigh-Jeans  formula),  we 
obtain  :•  :•  -  ;;  .  ••  ’•  " 


■"  .:V> '  12w2c2 

■•r  If  hw>  kT^j  then  S^T^.  ,  The  maximum  flow  of  energy  is 
obtained  in  the  case -  ©o  (air  'frequfencies  luminesce  [vysvechi- 
vayutsya])  then  S  ‘  s  <y  T4*  (o  =  5.67*  1CH>  erg/oa2 . s ec ...degree  is  the 
Stefan-Boltzmann  constant.  Besides  (l),  we  must  Still  fulfill  the 
condition  so  that  the  flow  of  heat  in  the  center  of  the  wave  is  >v 
equal  to  zero 


r2 


-  0 


2  - 


=  0. 


(2) 


We  shall  write  the  equation  for  heat  transfer  in  a  quiescent 

medium 

.11.  1  l  r?>HCTh  ;  x  „  clftl  .  . :  (3) 

r2  c£r  3 

Equation  (3)  is  equivalent  to' an  infinite  number  of  relation¬ 
ships  which  may  be  obtained  by  integrating  this  equation,  previously 
multiplied  by  r®  (m  2,  3>  k>-<>)-  Applying  the  method  of  moments, 
we  shall  satisfy  equation  (3)V having  required  approximate  fulfill¬ 
ment  of  only  two  integral  relationships,  corresponding  to  m  a  2  and  3 
Employing  (l)  and  (3),  we  obtain  the  equalities 

■  ~  i’4>  '  .  '  *  "  • 

l*Ka  ^  Tr2dr  =  Qo  -  J  SCT^)  r2  dt,  (h) 

0  0  '  ^  r 

r  l 

-1-  a  (  r^Tdr  =  -  r?  S(TA  )  -  cnr2  tQO  -f  2c0  (  rT'~  dr.  (5) 

dt  J  cb  T.  u  9.  4>  ,.?■  J  9 

0  0 

Despite  the  fact  that  out  of  an  infinite  number. we  have  satis¬ 
fied  only  two  equalities  (h)  and  (5),  we  can  hope  for  good  agreement 
between  approximate  and  accurate  solutions,  since  the  temperature 
behind  the  wave  front  changes  evenly  and  monotonlcally. 

Since  the  law  governing  the  distribution  of  temperature  be¬ 
hind  the  wave  front  is  basically  determined  fcy  the"  thermal  conduc¬ 
tivity,  we  shall  search  for  a  solution  in  a  form  similar  to  the 
progressing  problem  .....  : ^  , 

1  . .  'J"  ,<jS  4 ’ : • 

T  =  T0(t)  j*l  -  2,  j  k~J  (6) 

where,  however,  T  (t),  A(t),  r<k(t)  are  previously  unknown  functions 
satisfying  equalities  (l),  (U)  and  ($),  Applying  (6)  and  combining 
(2),  (h)  and  (5),  we  arrive  at  a  system  of  equations 
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A  ,  S  .  T<b  '  ,  r<fc  A  » 

r  -  / - =  0,  /  3  (  - - 7  ),  ■  0, 

A  ,  bT  ^  %  r  ;//;  *£  A  -  - 


•  (?) 


T 


2cq1<A2 


a(k  -  i)r2  0  . 3To 

<P.  : 


tJ"1  /  1°-  a  0,  a  T0(l  -  A*)*"1 


4. 


c  -  *  >■  •  ■■ '  ■  j ■  ■  \  .  ^  '  '  *  *"  !  ' 

Equation  (7)  is  easily  integrated*  In  the  Case  S  s 
Tire  obtain"  1  •'  "  • 


3k-l  „  ,  /n  mb-1  -  a  1 1<;Cl  x3(l-k)t  x  _  r«j>. 
:  =  Git  /  C2,  Tq  2kc0(3k  -  1)  “  T“ 


B-jc" 


3k-3n-2  ,  f  (i  -  a2)1"*  dA  /  Cv 


(8) 


Hie  constants  C^j  C2  and  are  determined  from  initial  conditions. 

If  initially  the  process  differs  from  that  described  by  ns,  but 
at  the  moment  of  time  tQ  there  is  formed  a  heat  waye  with  radius  rp 

and  a  distribution  of  temperature'  '  '  :'r / 

••  V-  !  1^»- 

21  k-1 


«*.’>■*  hstt 


then  the  initial  conditions  for  determining  the  three  constants 
T00>  r0>  Aqj  will  be 


T0(r0,  tQ)  =  T00,  r^  (t0)  =  rQ,  A(t0)  =  Aq. 

If  liberation  of  energy  occurs  at  a  point,  then  during  the 
initial  moment  solution  must  become  progressing.  Therefore,  we  must 
assume  :  V  :  . -  •'  .  '  ' 


co  Qo  k-1 

G2  =  c3  s  o,  Cl=  “  [  ITJlkT  3  » 


X (k)  * [  2k(3k-i)  3  f1  <  I  >f(  e~i  )  tr(  |  / ^ )] 1. 


(9) 


Let,  for  example,  n  =  0,  then  from  formulas  (8)  and  (9)  we 
obtain  the  series  of  expressions 

1  3k-l 

.  -  3k"2  „  -  r  °i  /3k-2 ,  3k~2 

A*  “  3k-l*  *<f>*  “  TWZTX  3  (3k-l3  > 


r 


rx=  A1  s  0, 


ri  =  -A  ' 


3 (k-1) 

a2  3  3k“l  ' 


(10) 


T  .  s 


a(k-l) 


-  T  -  f  >  '  i  r 

01  ~  4».l  "  1  2kOk--l;co  '  1 


°0 


X (3k-l) 
3  k  2 


T4>.max 

^.l 


3k-»l  1 

p  Bk^Hk-iT  ,  ¥-i 

(  3^-1  )  (3k  **  2)  ,  X 


Here,  the  values  corresponding  to  the  maximum  radius  of  the  front  of 
the  heat  wave  are  marked  by  an  asterisk,  and  the  values  at  the  moment 
when  all  energy  is  radiated  —  by  the  index  "l"}  T<$>  ,max  is  the 

maximum  temperature  at  the  wave  front,  A^  is  the  value  of  A  corres¬ 
ponding  to  this  temperature.  The  formulas  (10  )  describe  the  quite 
obvious  picture  of  heat  wave  propagation  with  consideration  for 
luminescence  in  the  motionless  gas. 

At  first  the  radius  of  the  wave  front  increases  in  time, 
reaching  its  maximum  value  rcjj,  after  which  the  dimensions  of  the 
heated  region  begin  to  shrink  and  the  radius  of  the  front  again 
decreases  to  zero,  at  which  time  all  the  original  energy  is  radiated 
"  [vysvechivayetsya].  The  temperature  at  the  wave  front  at  first  in¬ 
creases  from  zero  to  maximum  T<£  ma3c»  this  is  reached,  however,  after 
the  radius  of  the  front  begins  to  decrease.  .J  - 


The  extreme  values  of  all.  Quantities  are  easily  found  from 
the  formulas  (8)  as  well  as  for  any  arbitrary  value  of  n,  fnus,  for 
example*  for  A#  and  we  have 

.  i  .  r  ..  ^  n  .  . 

•  -  ‘  ..  .  .A*;.  ;  !  (  ;  i  ;  ;  ■  ■jzg  ;•  ; 

(3k  -  3N  -  2)  J  Z(A)dA  /  Z(A*)A#  a  0,  Z  =  (l  -  A2)  , 

1 

V  •.  ..  :  (11) 

i 

f  .  .  ,  ,  3(1  -  k)  o 

2A2  J  Z(A)dA  /  (1  -  A  )  y—%~rr  “  °*  . 

Ao  '  '  •  "  '  - 


We  note  that  from  the  formulas  (8),  in  the  case  when  tne 
liberation  of  heat  occurs  at  a  point, ' it  follows  that  the  tempera¬ 
ture  in  the  center  changes  according  to  the  progressing  law,  ana 
that  the  temperature  profile  coincides  with  the  progressing  pro^i  , 
which  is  clipped  off  at  a  distance  from  the  center  equal  «o  r^  s 
qA  (r^  q  is  the  radius  of  the  wave  for  the  progressing  problem). 

Solution  of  equation  (3)  could  be  sjught  in  another  form 


t  -  JLff  I c3^1] 
_  2kOk  -  1)Cq  J 


rA  \2l«(t) 


determining  r <J>  (t) ,  A(t)  and  0(  (t)  similarly  as  above  from  equa¬ 
tions  (l),  (ii)  and  (5).  ' 

Despite  the  fact  that,  mathematically,  the  formulas  (8)  de¬ 
scribe  the  propagation  of  the  heat  wave  in  all  stages  right  up  to 
complete  radiative  dissipation  at  a  certain  distance  *roo.  the  origin 
of  the  coordinates  the  solution  already  loses  its  physical  meaning 
bv  virtue  of  the  emergent  motion  of  the  gas.  Since  the  heat  wav  _ 
slows  down  quite  abruptly,  we  can  therefore  determine  qualitatively 
SfuSTof of  oar  solution,  as '  indicated^  Kompan^ets 
($),  from  the  conditions  -for  equality  between  tne  velocity  of  the  , 
wave  front  and  the  velocity  of  sound  in  its  center.  Energy  transfer 
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proceeds  further  by  gas-dynamical  means*  since  the  elementary  physi¬ 
cal  requirements  are  fulfilled  so  that  the  small  disturbances  over¬ 
take  the  front  of  pronounced  discontinuity.  Let  the  velocity  of 
sound  be  given  by  the  formula 


'0*  a  b  Ti 


1/2 


3Q 


Ti.  * 


v 


Since  the  temperature  profile  approaches  a  plateau  owing  to 
the  intense  heat  transfer*  Tj,  is  close  to  Tq}  we  then  obtain  une 

formula  which  determines  A3  from  the  condition  C*  =  ^ 


A3Z 


1  ‘(3k  -in  -  2 Tf  '  1  3  3(2n-l)  ,3f^n-lT  “  ft 


(12) 


ZA 


^jwsn-Tf 


p  -  JL  L* 
'  '  >  ? 


2n-l 

2~ 

)  -  B 


$ 


It  is  interesting  to  note  that  the  values  for  A3  approach 
unity  even  for  the  limiting  case  when  /3  ~  0.  These  values  are  de¬ 
termined  from  the  conditions  when  the  expressions  contained  in 
brackets  in  equation  (12)  are  reduced  to  aero.  In  the  case  n  -  0, 
expression  (12)  becomes  simple*  and  the  maximum  value  for  A^*  cor¬ 
responding  to  /3  s  0,  is  determined  by  the  formula 


3max  s  A# 


3k  -  2 

mi  a  him 

3k  -  1  * 


(13) 


The  velocity  of  the  front  of  the  heat  wave  at  this  time  be¬ 
comes  equal  to  zero.  For  /3  0*  the  values  for  A.,*  calculated 

from  formula  (12),  lie  still  closer  to  unity.  For  Example*  at  n  s  0, 

k  =  6,  y3=  200*  A«1  -  2.5  *10"^. 

Making  use  of  the  proximity  of  A3  to  unity*  expression  (12) 
can  be  linearized*  after  which  we  obtain 


-3max 


=  1 


k  -  n  -  1 

k  -  n  -  1  /  (k-1)  (3k  -  3n  -  27* 


(lli) 
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Knowing  the  temperature  distribution  along  the  radius  at  each 
moment.  It  ia.  easy  to  determine  the  share  of  radiated  energy 


The  results  of  computations  by  this  formula  for  the  case  k  s  3 
and  k  a  6  are  shown  in  the  figure.  We  notice  'that  the  amount  of 
radiated  energy  at  a  given  distance  depends  on  the  law  of  forced 
luminescence  only  through  the  value  A»  Therefore,  at  k  2  6,  n  s  0, 

^  0.9j  at  k  s  3,  — ■ —  ^  0.95. 

Q  Q 

For  the  case  when 

n  2  0,  k  =  6,^  =  200,  ^  Cr:  1  -  5  •  10**. 

From  the  above  it  follows  that  the  loss  of  fnergy  through 
radiation  from  a  wave  front  up  to  the  moment  when  r<£> *  eg  is  small, 
even  if  radiation  is  very  intense.  Physically,  this  results  from 
the  fact  that  an  increase  in  radiation  retards  the  wave  front,  and 
the  motion  of  the  gas  is  found  to  be  substantial  at  distances  closer 
to  the  origin  of  coordinates. 

In  conclusion,  let  us  note  that  in  the  case  when  the  length 
of  the  free  travel  of  radiation  is  commensurate  with  the  dimensions 
of  the  heated  region,  an  approximate  solution  can  be  found  by  the 
method  of  moments. 

The  author  extends  his  sincere  thanks  to  A.  S.  Kompaneyets 
for  hi3  valuable  assistance. 


1 

k-1 

y2)  dy.  (15) 


Q  -  Q  _  2 


t  37  $  *.  a 
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